Seasonal variation of atmospheric leptons as a probe of charm 
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The intensity of TeV atmospheric muons and neutrinos depends on the temperature in the strato- 
sphere. We show that the energy-dependence in the 100 TeV range of the correlation with tempera- 
ture is sensitive to the fraction of muons and neutrinos from decay of charmed hadrons. We discuss 
the prospects for using the temperature effect as observed in gigaton neutrino detectors to measure 
the charm contribution. 
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Atmospheric muons and neutrinos are important for 
neutrino telescopes, both as background in the search 
for high-energy neutrinos of astrophysical origin and be- 
cause they serve to calibrate the detectors. In addition, 
the atmospheric lepton spectra are of interest because 
the large size of current and proposed neutrino telescopes 
will make it possible to extend the measurements of at- 
mospheric muons and neutrinos beyond the reach of pre- 
vious measurements. It may finally be possible to achieve 
a long-sought goal of cosmic-ray science by detecting the 
contribution to the atmospheric leptons from decay of 
charmed hadrons, the "prompt" neutrinos and muons. 
At energies much lower than ~ 100 TeV, the prompt 
component is hidden by the much more abundant con- 
tributions from decay of charged pions and kaons. A 
temperature dependent signal of the charm contribution 
would complement the standard approach of looking for 
an isotropic component of the muon flux with a harder 
spectrum. 

It has long been known that the intensity of high- 
energy muons depends on the temperature in the strato- 
sphere p|. The variation arises from the corresponding 
expansion and contraction of the atmosphere at high al- 
titude where most production of the muons occurs. A 
similar effect is also expected for neutrinos for the same 
reason (2j. The temperature dependence of TeV muons 
measured in MINOS Q has recently been used as a probe 
of sudden stratospheric warmings 3 • IceCube reported 
a correlation of muon intensity with stratospheric tem- 
perature with high precision made possible by the high 
statistics of the kilometer-scale detector [j| . Ref. @ de- 
scribes how the temperature dependence of muon inten- 
sity in the 1 to 10 TeV range is sensitive to the kaon to 
pion ratio in the atmosphere. In this paper we point out 
that kilometer-scale neutrino detectors may have suffi- 
cient statistics in the 100 TeV range to be able to con- 
strain the contribution of charm decay to prompt leptons. 

Currently operating neutrino telescopes are Baikal Q , 
Antares Q and IceCube @. AMANDA [l(J, the prede- 
cessor of IceCube at the South Pole, was turned off at 
the beginning of 2009. With a planned volume of one 
cubic kilometer, IceCube will be by far the largest oper- 
ating neutrino detector. It is already running with 90% 
of its sensors installed and will be completed in 2011. 
There are plans for a still larger neutrino telescope in the 



Mediterranean Sea (Km3NeT) [Tl[, and there are also 
plans to expand the Baikal detector [l2T |. 

Preliminary measurements of muon [l3j and neu- 
trino [lij spectra with data from the partially con- 
structed (25% complete) IceCube in 2007 were re- 
ported at the International Cosmic Ray Conference, July 
2009. The atmospheric neutrino spectrum measured with 
AMANDA with data taken from 2000 to 2006 [h} al- 
ready extends to just above 30 TeV with high statistics. 
An unfolding analysis of the AMANDA [15j data extends 
the reach of the same data set to above 100 TeV. No con- 
tribution of a prompt component is yet visible in the un- 
folded neutrino spectrum. The full IceCube will be large 
enough to extend measurements of atmospheric neutri- 
nos well beyond 100 TeV with good statistics. From the 
point of view of precise analysis of the energy dependence 
of the correlation coefficient, the extremely high statis- 
tics available with muons in a kilometer scale detector is 
also important. 

Fluxes of secondary cosmic-ray leptons in the atmo- 
sphere can be described to a good approximation by a 
set of formulas in which each element corresponds to one 
of the processes involved in their production [ll| [TtJ • For 
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where 4>n{Eu) = dN / d\n{E v ) is the primary spectrum 
of nucleons (N) evaluated at the energy of the neutrino. 
The three terms in brackets correspond to production 
from leptonic and semi-leptonic decays of pions, kaons 
and charmed hadrons respectively. The equation for 
muons is similar at energies sufficiently high (> 100 GeV) 
so that energy loss and decay of muons in the atmosphere 
can be neglected. The numerator of each term is of the 
form 
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with i = 7T^~ , K, charm and BR{ V is the branching ratio 
for i — > v. The first Z-factor in the numerator is the 
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spectrum weighted moment of the cross section for a nu- 
cleon (N) to produce a secondary hadron i from a target 
nucleus in the atmosphere, and the second Z-factor is 
the corresponding moment of the decay distribution for 
i — > v + X. Thus, for example for N + air ->• K + + X, 
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where x = E K + /En and 7 » 1.7 is the integral spectral 
index of the incident spectrum of cosmic-ray nucleons. 
The denominator of Eq. [2] is the ratio of nucleon inter- 
action to attenuation length. For calculations below we 
use numerical values of the parameters from Ref. [l6| . 
Parameters for charm production are discussed below. 

The denominator of each term in Eq. Q] reflects an 
important physical property of meson decay in the 
atmosphere-the competition between decay and reinter- 
action of hadrons. The critical energy depends on the 
zenith angle 9 of the cascade in the atmosphere, and is 
of the form 
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where 9* is the local zenith angle at lepton produc- 
tion taking account of the curvature of the Earth [l7| . 
The Ci are characteristic energies for each channel and 
ho ~ 6.4 km is the scale height in an exponential approx- 
imation to the density of the atmosphere at high alti- 
tude. When Ei < ej/ cos 9 the mesons decay so that the 
low-energy neutrino spectrum has the same power law 
index as the primary cosmic-ray spectrum. At high en- 
ergy, the contribution of each term gradually steepens so 
that asymtotically, the neutrino spectrum is one power 
steeper than the primary spectrum. The characteristic 
energies are 

e w = U5 GeV, e K ± = 850 GeV, and 

Ccharm ~ 5 X 10 7 GeV. 

These differences lead to a characteristic pattern of 
the contributions of the various channels to the flux of 
neutrinos: first the pion contribution steepens, then the 
kaon contribution and finally (at much higher energy) the 
contribution from decay of charmed hadrons becomes the 
main source of neutrinos and muons in the atmosphere. 
In addition, the suppression of the charged pion contribu- 
tion to the neutrino flux starts early because the quantity 
B-nv = 2.8 in the denominator of the first term of Eq. [T]is 
anomalously large compared to the corresponding factors 
for strange and charmed hadrons. This is a consequence 
of the kinematics of tt^ — s> ^ + v decay in flight in which 
the muon carries most of the energy because its mass is 
comparable to that of the parent pion. 

Fig. [T] shows the fractional contribution of the main 
hadronic channels to the production of neutrinos and 
muons in the atmosphere. Kaons actually become the 
dominant source of neutrinos above ~ 100 GeV. For 
muons, the dominant channel is never kaons, but the 
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FIG. 1: Top: fraction of muon- neutrinos from pion decay, 
kaon decay and charm decay (RQPM [241 ] ) as a function of 
neutrino energy; Bottom: same for muons. Solid lines are for 
vertical and dashed lines for 60° (see text for discussion). 



contribution of kaons does increase significantly in the 
TeV energy range. An interesting manifestation of the 
increasing importance of kaons at high energy is the mea- 
surement of the charge ratio of atmospheric muons with 
energies at production in the TeV range 0, The 
ratio is observed to increase from its value of 
1.27 around 100 GeV to w 1.37 above a TeV. The in- 
crease is attributed to the importance of the associated 
production of kaons, 
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which makes the +/— charge ratio larger for kaons than 
for pions. Production of strange particles and anti- 
particles-in particular the process in Eq. [5]-is highly 
asymmetric in the forward fragmentation region. Pro- 
duction of A is favored because it has constituents in 
common with the valence quarks of the proton. 

At some energy, perhaps around 100 TeV or somewhat 
above, the decay of charmed hadrons will become the 
main source of all atmospheric leptons. Where the tran- 
sition occurs is uncertain because of the wide variation in 
the literature on the level of charm production at large 
Feynman x. Like associated production of strangeness 
via Eq. [5J the production of is also highly asymmetric 
in the forward fragmentation region. A measurement by 
the Fermilab Experiment 781 (SELEX) [lc| shows that 
(cr c - at-) I (<7 C + o-c) ~ 1 for 0.2 < x < 0.7 for 600 GeV 
protons on a fixed target. Here a c and a s represent re- 
spectively the production of A+, which does, and AjT, 
which does not have valence quarks in common with the 
beam proton. A large contribution from the charmed 
analog process to Eq. [5] could be classified as " intrinsic" 
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charm 21]. There is some support for a component of 
intrinsic charm from recent measurements of charm pro- 
duction on different nuclear targets [12], as discussed in 
Ref. H. 

The RQPM model of charm production [24[ includes 
such a contribution and predicts a relatively high level of 
prompt leptons, as illustrated in Fig. Q] Calculations us- 
ing perturbative QCD tend to give lower levels of charm 
production. As an example, a recent calculation within a 
perturbative QCD framework (25j predicts a contribution 
from charm decay roughly an order of magnitude lower 
than the RQPM model. In what follows we evaluate the 
sensitivity of seasonal effects to prompt atmospheric lep- 
tons for three different assumptions about the level of 
charm production. We compare RQPM with charm pro- 
duction at the level of Ref. [25[ (henceforth ERS) and 
with a somewhat arbitrary intermediate model described 
in terms of the parameters of Eq. [5J Representing the 
weighted sums of the various charm channels and their 
scmi-leptonic decay modes with i?N-charm = 5 x 10~ 4 , 
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= 0.13 and an effective semi-leptonic 
branching fraction for charm decays of 0.14, gives a level 
of charm production about a factor of two lower than 
RQPM. 




FIG. 2: Differential correlation coefficient as a function of 
lepton energy for near vertical muons and muon neutrinos. 
Solid lines: no charm; dashed lines: with RQPM charm (24[. 



In an isothermal approximation, the density of the at- 
mosphere is described by an exponential with a scale 
height of ho « 6.4 km, where the numerical value is appli- 
cable to the stratosphere where most high-energy muons 
and neutrinos originate. From the ideal gas equation re- 
lating density and pressure, one finds ho — RT. With 
this relation it is then possible from Eqs.|3]and[T]to cal- 
culate the variation of the flux with temperature. The 
deviation of the atmosphere from isothermal is accounted 
for to first order in 5T/T by weighting the actual tem- 
perature profile by the profile of production of the atmo- 
spheric leptons. It is conventional to define a correlation 
coefficient, which relates the fractional change in the 
atmospheric muon intensity to the fractional change in 
the temperature. One can describe the atmospheric neu- 
trinos in parallel, with a coefficient a v [2j . In differential 
form the correlation coefficient is 
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The derivative can be calulated directly from the expres- 
sions [T] and IH with the results shown in Fig. [2l 

The features of the curves in Fig. [2] correspond to the 
properties of the various terms in Eq. [TJ At extremely 
high energy, the second term in the denominator dom- 
inates, which means the the intensity of muons is pro- 
portional to temperature and the correlation coefficient 
is unity. Asymptopia is reached at lower energy for pi- 
ons than for kaons, so the inclusion of kaons delays the 
increase of the correlation coefficient so somewhat higher 
energy. This is more apparent for neutrinos than for 
muons because kaons give a larger fractional contribu- 
tion to neutrinos than to muons. The increase of the 
correlation coefficient is delayed to higher energy at large 
zenith angle because the critical energy (Eq. 2J is larger. 
Since the critical energy for the charm contribution is at 
extremely high energy (~ 5 x 10 7 GeV), the prompt lep- 
tons from charm decay have no temperature dependence 
until ~ 10 7 GeV. Thus, a significant contribution from 
charm will suppress the correlation coefficient above an 
energy that depends on the magnitude of the charm con- 
tribution to lepton production. For the example shown 



in Fig. [5]in which the RQPM model of Bugaev et al. [24 1 
is used, the suppression begins already between 10 and 
100 TeV depending on zenith angle. 

Measured rates depend on the convolution of the lep- 
ton spectrum with the detector response and effective 
area, which depends on lepton energy and direction. To 
compare with measurements, it is therefore necessary to 
make the following calculation: 
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Here the subscript i = u, /i and Ai >e s (E, 9) is the energy- 
dependent effective area of the detector as seen from a 
direction 9. In the remainder of the paper we explore 
some examples relevant for IceCube. 

To make an estimate for muons we assume a step func- 
tion at energies corresponding to the minimum needed 
to penetrate through IceCube treated as a sphere with 
a projected area of one square kilometer. Table Q] gives 
some rates for muons with zenith angle < 30°. The cor- 
relation coefficient is 10% lower with RQPM charm for 
Efj, > 10 TeV than with no charm and the rate is more 
than one Hz per square kilometer. For E^ > 100 TeV 
the correlation coefficient is only 0.5 with charm while 
the event rate is much lower, but still of order 10 5 events 
per year in a kilometer scale detector. 

For neutrinos we use the effective areas computed for 
IceCube-40 [26| multiplied by a factor of two to estimate 
the event numbers for a full cubic kilometer detector. 
Table |H] gives correlation coefficient and expected rates 
for three different ranges of zenith angle for neutrino- 
induced muons coming into the detector from below. The 
three ranges of zenith angle are 1) 90° < 9 < 120°, 2) 
120° < 9 < 150°, and 3) 150° < 9 < 180°, which have 
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0.89 
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0.0025 


0.53 


0.0049 


0.91 


0.0028 


0.71 


0.0036 



TABLE I: Correlation coefficients for muons with (8 < 30° 
for three levels of charm (energy in TeV; rate in Hz/km 2 ). 



£,, min (TeV) 


no charm 


RQPM charm 


Zone 1 
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Events /yr 
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Events /yr 


all 


0.54 


16000 


0.52 


17000 


3 


0.70 


5900 


0.62 


6300 


30 


0.94 


350 


0.72 


450 


£,, min (TeV) 


no charm 


RQPM charm 


Zone 2 


a 


Events/yr 


a 


Events/yr 


all 


0.66 


6000 


0.62 


6400 


3 


0.88 


1230 


0.75 


1450 


30 


0.98 


37 


0.46 


80 


£,,min(TeV) 


no charm 


RQPM charm 


Zone 3 


a 


Events/yr 


a 


Events/yr 


all 


0.68 


1650 


0.64 


1750 


3 


0.91 


260 


0.75 


320 


30 


0.99 


5.2 


0.41 
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TABLE II: Correlation coefficients with and without charm 
for neutrinos in three zones of the atmosphere (see text). 



solid angles respectively of 7T, 0.737T and 0.277T. For a 
detector at the South Pole, zone 1 corresponds approxi- 
mately to the Southern temperate atmosphere, zone 2 to 
the tropics and zone 3 to the Northern temperate atmo- 
sphere. Absorption of neutrinos in the Earth is significant 
for the high-energy entries in the vertical bins. 



For lepton energies below 10 PeV the contribution from 
charm decay does not depend on temperature because 
charmed hadrons decay before interacting as a conse- 
qence of their short decay times. In contrast, by 100 
TeV pions and kaons are fully asmyptotic in the sense 
that their decay probability is proportional to tempera- 
ture. This situation leads to the possibility of inferring 
the level of prompt lepton production from a precise mea- 
surement of the correlation of the rate of leptons with 
stratospheric temperature in the 10 to 100 TeV energy 
range. There are several practical problems that must be 
overcome to realize this possibility. First, a high-energy 
event sample must be defined with sufficient energy res- 
olution and sufficient statistics to see the effect of the 
charm contribution on the correlation with temperature. 
In the case of muons, there will be sufficient statistics 
in a kilometer scale detector, but the energy assignment 
of individual muons may be complicated by the presence 
of accompanying muons from the same cosmic-ray pri- 
mary. In the case of neutrinos it will be necessary to in- 
tegrate over a large fraction of the surface of the Earth, 
although about 25% of up-ward atmospheric neutrinos 
are generated in the atmosphere above Antarctica, where 
the stratospheric temperature variation is strongest. The 
practical application of this probe of charm is being de- 
veloped in association with analysis of measured tem- 
perature effects in IceCube and is most promising for 
high-energy muons because of the large statistics. 
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